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12/04/2011 ABSTRACT 



Analysing all Galaxy and Mass Assembly (GAMA) galaxies within a factor two 
(±0.3 dex) of the stellar mass of the Milky Way (MW), there is a 11.9% chance that 
one of these galaxies will have a close companion (within a projected separation of 
70 kpc and radial separation of 400 km/s) that is at least as massive as the Large 
Magellanic Cloud (LMC). Two close companions at least as massive as the Small 
Magellanic Cloud (SMC) are rare at the 3.4% level. Two full analogues to the MW- 
LMC-SMC system were found in GAMA (all galaxies late- type and star forming), 
suggesting such a combination of close together, late-type, star-forming galaxies is 
rare: only 0.4% of MW mass galaxies (in the range where we could observe both the 
LMC and SMC) have such a system. In summary, the MW-LMC-SMC system is a 
2.7a event (when recast into Gaussian statistics). 

Using cross-correlation comparisons we find that there is a preference for SMC- 
LMC binary pair analogues to be located within 2 Mpc of a range of different lumi- 
nosity groups. There is a particular preference is for such binaries to be located near 
LG luminosity systems. When these groups are subdivided into small magnitude gap 
and large magnitude gap subsets, the binaries prefer to be spatially associated with 
the small magnitude gap systems. These systems will be dynamically less evolved, but 
still offer the same amount of gravitational dark matter. This suggests that binaries 
such as the SMC-LMC might be transient systems, usually destroyed during vigorous 
merger events. Details of a particularly striking analogue to the MW-SMC-LMC and 
M31 complex are included. 
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1 INTRODUCTION 

The Local Group (LG), and more specifically the Milky- Way 
(MW) is the most thoro ughly explored dark matter com- 
plex in the Universe (e.g. lMateolll998l ; Ivan den BergblliocT 



se (e. : . . 

iBenson et al.1 |2002| ; iKarachentsev et al.1 120091 : iFont et al 
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l201ll) . However, question marks remain over how typ- 
ical the MW halo is in the context of the Uni- 
verse and how unusual its ga l axy occupation statistics 
are (e.g. [Bovlan-Kolchin et al.1 l201ll ; lLovell et all l201ll : 
iTollerud et al.ll201ll ; IWeisz et alj|201 lK . It is important that 
we fully understand how representative the MW halo is 
since, by virtue of proximity, it will always be the en- 
vironment that will contain the faintest known galaxies 
and the broadest range of galaxy masses. It will also be 
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the halo from which we can derive the most informa- 
tion about its formation history. Knowing which satellites 
of the MW halo are typical within similar-mass similar- 
redshift haloes will either severely tighten or relax the pre- 
dictive requirements of N-body semi-analytic galaxy for- 
mation codes. Currently it is acknowledged that simula- 
tions struggle to predict the full distribution of MW satel- 
lite galaxies, these problems are particularly manifest for 
the brightest satellites: Large and Small Magellanic Clouds 
(LMC and SMC) (e . g. IBenson et all |2002| ; IKoposov et~all 
l2009l ; IOkamoto et al.ll2010l ). 

We are set to learn a vast amount about the MW 
halo in the coming d ecades. In the near future GAIA 
l|Wilkinson et all I2005T ) will measure space motions and 
properties for 2 billion stars in the LG which includes all 
known member galaxies. Amongst likely discoveries, we will 
learn about dynamical equilibrium, or lack of it, for the first 
time. Building up to these hugely detailed surveys it is im- 
portant we discover where the MW halo fits into the bigger 
picture. Only then can we apply what we know about the 
MW to the ACDM (or some variant) model of the Universe. 
Combining near-field cosmology (LG scale) and far-field 
cosmology (redshift surveys) is key to completing the full 
pictu re of galaxy formation (|Freeman Sz Bland-Hawthornl 
12002ft . 

This work puts the investigation of the MW halo into 
an observational cosmological context by using data from 
the Galaxy and Mass Assembly project (GAM A). GAMA is 
a multi-wavelength photometric and spectroscopic survey, 
and was designed to answer questions about how matter 
has assembled on a huge variety of scales: filaments, clus- 
ters, groups and galaxies. The first phase of the redshift 
survey was conducted on the AA T (known as GAM A-I) and 
these data are used in this work (|Driver et alll201lT ). In this 
work we use GAMA redshifts to search for close compan- 
ions to MW mass galaxies. These systems will be Milky- 
Way Magellanic Cloud Analogues (MMAs from here). We 
use this sample to construct statistics on the rarity of SMC 
and LMC type (star forming late-type galaxies) close com- 
panions around L* late-type moderate star-formation rate 
spiral galaxies like our own MW. 

In Section [2] we discuss the data used in this work in 
detail. In Section [3] we present the statistics for finding 
MW-LMC and MW-LMC-SMC type systems, allowing us 
to quantify the apparent rarity of MW like systems. In Sec- 
tion U we investigate the environment that SMC-LMC type 
binaries are most commonly located in, and relate this to 
some of the defining characteristics of the LG. 

Data for the LG was calculated using distance indica- 
tors without any Ho dependence. As such it is appropriate 
(and consistent with the main body of LG literature) to 
convert GAMA data into true distance. To make the appro- 
priate conversions we take the latest WMA P 7 cosmology: 
p M = 0.27, Q. A = 0.73 and H = 70kms _1 (jKomatsu et all 

2 DATA 

The Galaxy and Mass Assembly project (GAMA) is a 
major new multi- wavelength spectroscopic galaxy survey 
l|Driver et al.ll201lft . The first 3 years of data obtained is re- 
ferred to as GAMA-I, and is the data used for this work. 
The GAMA-I data used here contain ~86K redshifts to 
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Table 1. Data used in this work. Data origin: Ivan den Bereli 
l|2000h ". assumed M31 values for MW property 6 . | Flynn et al. 
[|2006h c , Nasa Extr agalactic Database**, I James fc Ivorvl <201llv C 
iNichols et al ] I^Olli y . All magnitudes use are intrinsic, i.e. Hq = 
70kms _1 . M r is derived using My and the B — V colour using 
the relevant Lupton SDSS photometric transform equation . 



r — 19 A over ~144 deg 2 , with a survey design aimed 
at providing an exceptionally unifo rm spatial completeness 
and high close pa i r completeness llRobotham et al . 2010a; 
iBaldrv et al]|20ld : iDrlver et alll201lft . 

Extens ive details of the G AMA survey characteristics 
are given in lDriver et al. |201lft , with the survey input cata- 
logue described in Baldrv et all ||2010|) a nd the spectroscopic 
tiling algorithm in iRobotham et al. I d2010ah. Additional 
data used for this w ork are stellar masses ILTavlor et al.ll201lh 
and G 3 Cvl groups l|Robotham et al.ll201ll l 

Table [T] shows the important MW, LMC and SMC 
values used for this work including distance and ra- 
dial/tangential velocity separations to the MW. Conversions 
to the native GAMA r-band are made using the Lupton 
transforms Q 

3 FINDING MW MAGELLANIC CLOUD 
ANALOGUES 



There are a number of questions that could be stated with 
respect to investigating MMAs. For clarity this Section an- 
swers the following: 'how common is it to observe 1 or 2 
relatively large (stellar mass > 1 x 1O 8 M0), star-forming 
satellites close to galaxies with stellar masses within 0.3 dex 
of the MW?' While this question lacks simplicity, it is at 
least possible to answer it in a meaningful, and reproducible, 
manner. This is a pertinent question to ask, since galaxy for- 
mation models have trouble replicating the presence of very 
bright satellit es in close proximity to simulated galaxies like 
the MW (e.g. IBenson et alj|2002» . 

First we define 'similar' stellar mass to me an within 0.3 
dex o f the MW mass (M s -mw = 5 x 10 10 M iFlvnn et all 
2006). Similarly we have to quantify 'close'. A number of 
studies have been contacted that allow us to estimate the 
SMC/LMC distance and radial/tangential velocity compo- 
nents. Much endeavour has been invest ed in proper mo- 



tion m e asurements by multiple teams dKallivavali l et al 



2006al lbT; ICosta et all 120091 : IVieira et all 120101 : ICosta et al 
20 111 , e.g.) which has allowed the more difficult to ob 



serve tangential velocity components to be estimated. These 
numbers have been further refined by detailed simulations 
citep[e.g.]bekk08,nichll. The SMC is ~ 60 kpc from the MW 
and travelling radially away at ~ 17kms _1 and 301 kms -1 
tangentially with respect to the MW: a net velocity of 
302 km s" 1 l|Nichols et alll201lft . The LMC is - 50 kpc from 
the MW and travelling away at ~ 89kms _1 radially and 
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367 km s tangentia lly with respect to the MW: a net ve- 
locity of 378 km s" 1 (|Nichols et alj|201lh . To conservatively 
recover all systems where the galaxies are in such spatial- 
velocity proximity, we create a catalogue of pairs for this 
work where the projected separation is r aop _p ro j < 70kpc 
and the radially velocity separation is i> SC p_ ra d < 400 km s _1 . 
It should be stressed that whilst varying the precise defini- 
tions of 'close' and 'similar' does impact the total number of 
systems recovered and the fractions, it has minimal impact 
on the ratios between different population fractions. 

A consequence of such a selection, and any similar, 
is that we do not truly distinguish between systems that 
have close pairs where all three galaxies show signs of inde- 
penden t 3-body interac tions (e.g. the M81-M82-NGC3077 
group; lYun et al.l 119941 ) and systems like the MW -LMC- 
SMC that had a binary i nfall formation history (e.g. iBekkil 
120081 ; iNichols et ai]|201lf ). In data of GAM A quality there 
is too little phase space information (two dimensions of high 
accuracy spatial positions and one dimension of low accuracy 
velocities) to constrain the likely formation history of any 
system given the selection criteria stated above. Therefore, 
any fractions quoted should be considered as upper limits for 
finding systems that had a similar binary infall history to 
the MW-LMC-SMC but real limits for finding systems with 
exceptionally small phase separations regardless of the spe- 
cific formation mechanism (e.g. both MW-LMC-SMC and 
M81-M82-NGC3077). 

A second point to consider is how sensitive we are 
to the instantaneous flux of the SMC and LMC, both of 
which have had compli cated star formation histories (e.g. 
lHarris fc ZaiitskvllioM ). Since the complicated tidal inter- 
actions between the MW-LMC-SMC are known to trigger a 
large amount of the star formation (e.g. IZaritskv fe Harris! 
2004), it seems prudent to use the current luminosities and 
phase positions of all galaxies concerned in order to find 
analogues. This simplifies the search compared to looking 
at larger distances because such systems might have entirely 
different stellar mass content due t o experiencing a more qui - 
escent evolutionary history. In fact lZaritskv fc Harris! (|2004l ) 
suggest that as much as 70% of the stellar content of the 
SMC may have been formed due to interaction triggered 
star formation. 

We use the r < 19.4 GAMA-I survey data. Applying 
these selection limits to recover all systems that have similar 
pairwise properties to the MW and the Magellanic Clouds 
creates a catalogue containing 3,731 galaxy-galaxy pairs and 
6,840 unique galaxies with no redshift limits applied. Obvi- 
ously some galaxies will have more than one pair (the MW 
has two — the SMC and LMC). We create complexes that 
contain all galaxy-galaxy associations and count this as a 
single 'pair' system, i.e. the MW-LMC-SMC as a single 'pair' 
system. Throughout we use the beta distri bution to put ro - 
bust estimates on the sampling statistics l|Cameronll201 ll ). 
giving us firm limits on the 68 percentile probability range 
where relevant. 

In the LMC depth sample (0.01 < z pair < 0.089 for 
r < 19.4 mag) we find 1,642 galaxies that have M s = 
Mb-mw within a 0.3 dex range. Of these 286 galaxies are 
the dominant galaxy in the system (240 pairs, 39 triplets, 
6 quartets, 1 quintet), where there are 340 minor galaxies 
and 626 galaxies in total. This suggests that 17.4% of MW 
mass galaxies have at least one 'close' companion. Of these 



paired systems 56/286 have late- type dominant pair galaxies 
(19.6%). Of all the systems that have a late-type brightest 
pair galaxy (BPG), 34 of the minor pair galaxies are late- 
type. Since complexes can contain more than one pair, this 
translates as 31/56 late-type BPGs have at least one late- 
type minor companion (55.4%) and 30/56 that have 100% 
late-type companions (53.6%). 

Of these 30 systems that are 100% late-type, all have 
some observable Ha emission in the larger galaxy (100%) 
and 32/34 minor companions have some amount of observ- 
able Ha emission (94.1%). Comparing to the larger sam- 
ple of MW mass selected galaxies this means that 30/1,642 
(1.8%, beta range 1.5%-2.2%) of all MW mass galaxies are 
late- type and in a pair, where all of the minor galaxies are 
late- type and all galaxies are star forming. Thus, approxi- 
mate analogues of the MW system (where we just require 
all close companions to be Magellanic Cloud like) are rare 
at the less than 2% level. In this cascade of fractions, the 
most unusual characteristic is to find a late-type MW mass 
galaxy in a pair at all, followed by the dominant galaxy be- 
ing late-type given that it is in a pair. Once these criteria 
have been met the chance of finding a late-type compan- 
ion, and star formation in both galaxies, is remarkably high 
(over 50%). This means the discussion of the uniqueness of 
the MW system is largely driven by how rare galaxy pairs 
are at low redshift. 

Three triple systems are present in the final selection. 
The effective stellar mass limit is less well defined than 
the r-band limit (which has to assume an approximate k- 
correction), but based on when the number counts begin to 
turn over for galaxies with stellar mass less than the SMC 
the survey is complete out to z ~ 0.055. The total observ- 
able volume, using the standard cosmology of this paper, is 
1.8 x 10 5 Mpc 3 . Two of the three triple systems fall within 
this redshift range. Since 414 MW stellar mass ±0.3 dex 
galaxies are within this redshift limit, full analogues of the 
MW-LMC-SMC system are rare at the 0.4% level (0.3% to 
1.1%). In terms of space density, we find 1.1 x 10" 5 Mpc" 3 
full analogues in GAMA (in a volume of 1.8 x 10 5 Mpc 3 ). It 
is 6 times more likely that a MW mass galaxy with two SMC 
mass companions will have early-type morphology. Casting 
this figure in terms of Gaussian statistics, full analogues 
of the MW halo are rare at the 2.7cr level when searching 
around L* galaxies. 

Of the two systems, only one has a minor companion 
close to the stellar mass of the SMC (the other has two LMC 
mass companions). This system is the nearest analogue to 
the MW system found in the GAMA database, possessing 
a dominant star forming spiral galaxy with A4 S = 3.1 x 
10 10 M o (Ms-mw = 5 x 1O 1O M ), a more massive late- 
type companion with M B — 6.1 x 10 9 Mq {M s ~lmc = 
2.3 x 10 9 Mq) and a less massive late-type companion with 
Ms = 6.1xlO 8 M {Ms-smc = 5.3 x 10 s M ). Both of the 
companions are more massive than their Magellanic Cloud 
equivalents, but the smaller one is very close to the mass of 
the SMC. 

While the numbers found in this work suggest the MW- 
Magellanic Cloud system is cosmologically rare, we at least 
know that such a combination of galaxy diagnostics is not 
entirely unique. It is interesting to note that this system 
does not find the two companions to be in a close binary 
(as the LMC and SMC appear to be). This is similar to the 
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r/Mpc r/Mpc 

Figure 1. Left panel shows the cross-correlation between close binaries that are analogues to the SMC / LMC pair (Magellanic Cloud 
Analogues: MCAs) and G 3 Cvl groups of different extrapolated flux content within 0.5 mag of the value stated in the legend. The total 
r-band magnitude of the Local Group is M t —lq =~ —22.2, and the total magnitude of the MW halo is M, — mw =~ —21.3. The right 
panel splits the most highly clustered MCA-group subset (SM r = —22.5) into two further subsets: one with M r _ 2 nd — M r _ 1B t less than 
the median magnitude gap value of 0.7 for the groups selects, and another with M r _ 2 nd — M r _ lB t larger than the median. In both plots 
the y-axis is the standard cross-correlation excess, and is the relative excess compared to the random volume. 



findings of I James fc Ivorvl f|201 lh who, using an Ha limited 
imaging survey of 143 spiral galaxies, did not find a single 
MW-Magellanic Cloud analogue that had two companions 
in a close binary formation. This suggests that the binary na- 
ture of the Magellanic clouds might be their defining unique 
feature. 

To c ompare to the analysi s of th e Millenium II simu- 
lation bv iBovlan-Kolchin et ail l|201ll ). we select the 1,642 
galaxies that are within 0.3 dex of M s -mw as discussed 
above. We now apply a mass selection on the pairs. If we 
state that any minor pair galaxy must be the mass of the 
LMC or larger (M s -lmc = 2.3 x 10 9 M ) then we find 196 
galaxies that have a sub-halo occupied by a close companion 
that is at least that massive. This means that given a halo 
that is occupied by a MW mass galaxy, there is a 11.9% 
(11.2%-12.8%) chance that a galaxy at least as massive as 
the LMC occupies a sub-halo. 

IBovlan-Kolchin et al.l (|201lf ) find that in haloes that 
have a MW mass galaxy there will be a sub-halo containing 
a galaxy at least as massive as the LMC 3-8% of the time if 
the halo is ~ 1 x 10 12 M and 20-27% of the time if the halo 
is ~ 2.5 x 10 12 Mq. The median G 3 Cvl halo mass we find 
for MW mass galaxies is ~ 2 x 10 12 Mq (using the functional 
A scaling in Robotham et al. 2011). This num ber is consis- 
tent with the halo mass of the MW given bv iLi fc White! 
(2008). Our probability range of LMC or more massive sub- 
haloes (11.2%-12.8%) is als o between the ranges stated by 
IBovlan-Kolchin et al.l (|201lh . suggesting it is broadly con- 
sistent with their results. 

To calculate similar statistics for high order systems we 
apply the 0.01 < z < 0.055 SMC stellar mass depth limit 
(leaving 414 galaxies), and remove galaxies with stellar mass 



less than the SMC. We find 14/414 MW mass galaxies that 
have 3 or more galaxies in the pair system, i.e. 3.4% (2.7%- 
4.5%) of systems where we expect to be able to observe 
both the LMC and SMC have at least 2 galaxies with stellar 
mass greater than the SMC. The constraint is less tight due 
to poorer number statistics, a consequence of the smaller 
parent population within the SMC observable redshift limit. 
Purely framing the discussion in terms of how likely it is 
for baryons to occupy sub-haloes, these results imply that, 
assuming the dominant galaxy in the halo is similar in mass 
to the MW, it is 3.5 times less likely to have two sub-haloes 
with at least SMC stellar mass compared to 1 sub-halo with 
at least LMC stellar mass. 

4 WHERE DO MAGELLANIC CLOUDS 
ANALOGUES LIVE? 

An actively discussed mechanism for explaining the pres- 
ence of the Magellanic Clouds in the MW halo is binary 

aL||2009; 



infall (see discussion in 
lYang fc Hammeijl2010l ; 



Bekkill2008l: 



Nichols et al 



iallivavalil et _.. _ 
20 111 ). This model 



sumes the SMC and LMC were a binary pair that entered 
the MW-halo / LG-complex simultaneously. Observation- 
ally we can determine the viability of such a mechanism by 
searching first for close binary analogues to the Magellanic 
Clouds, and then determining the cross correlation these 
systems have with various luminosity grou ps taken from the 
GAM A Galaxy Group Catalogue (G 3 Cvl. iRobotham et all 

To select close binary Magellanic Cloud Analogues 
(MCAs) we search for all galaxies that have a close com- 
panion within a projected separation of 100 kpc and veloc- 
ity separation of 100 km/s, where both galaxies are between 
— 19 < M r < —17. This selection conservatively selects all 
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binary pairs with major characteristics in common with the 
SMC and LMC. 46 such binary MCA systems are found, 
where 1,929 galaxies fall inside the redshift and magnitude 
selection limits. This implies 4.8% (4.3%-5.3%) of galaxies 
in the specified magnitude range are in such close binary 
systems. 

Taking the G 3 C vl catalogue we use the cross- 
correlation approach of lCroft et alj (| 1999T ) to determine how 
spatially associated the MCAs are with different luminosity 
groups (within 0.5 mag of the stated values), where we use 
the extrapolated r-band flux content given in the G 3 Cvl 
catalogue. Errors are estimated through jack-knife resam- 
pling (measuring the variance in the cross correlation signal 
when binary pairs are excluded) and random volume cones 
are generated through uniform pointing (within the redshift 
range explored the sample is complete) . The left-hand panel 
of Figure [T] shows how closely associated the MCAs are to 
different types of groups. The medium galaxy flux plotted on 
the left panel (green line) includes groups with the same flux 
content as the entire LG. The red line shows more massive 
systems than the LG, the blue line includes groups with 
the same total flux as just the MW halo. It is clear that 
MCAs are more likely to be found in proximity to LG mass 
complexes, suggesting that the presence of M31 near to the 
MW halo has enhanced the probability of observing a bi- 
nary system like the SMC-LMC. A caveat to this result is 
the cross-correlations are affected by the subset comparisons 
chosen, and given the known interplay between luminosity 
and interactions the systems detected at large separations 
are likely to be in a different evolutionary state to those at 
close separations. 

To further investigate the effect seen for EM r = —22.5 
groups this cross-correlation was subdivided into two sub- 
sets: one where the magnitude gap of the brightest two 
galaxies (M r _ 2 ,id — M r _ ls t) is larger than the median value 
of 0.7 (so the central galaxy dominates) and another where 
the magnitude gap is smaller than the median (so the central 
galaxy does not dominate). The right-hand panel of Figure [T] 
shows the cross-correlations for these two subdivisions, with 
the original full sample plotted also. There is a marginally 
significant preference for MCAs to be more spatially associ- 
ated with groups with small magnitude gaps. These systems 
should be dynamically less evolved since the central galaxy 
is not as dominant within the group, and in fact is one of 
the major characteristics of the LG: the MW and M31 are 
similarly massive galaxies that would be found in the 'small 
gap' sample. This data lends evidence that the MW halo 
should not be considered in isolation when determining the 
occupation probability of the Magellanic Clouds. The poten- 
tial role of M31 on the presence of the LMC and SMC near 
the MW has been considered in recent simulatio n work (e.g. 
iKallivavalil et aL I l2009l ; lYang fc Hammerll2010l ), suggesting 
that its presence might be significant rather than coinciden- 
tal. 

To quantify this effect differently we move from con- 
sidering the cross correlation of MCAs with groups, to di- 
rectly searching for L* pairs in close proximity. We create 
the L* pair sample by searching for all galaxies between 
-21.9 < M r < -20.9 (the MW and M31 fall well inside this 
selection), within r pro j < l,000kpc and v Bep < 500 km s" 1 . 
Of the 302 galaxies that fall within the magnitude selection 
96 (32%) have a close V companion (i.e. there are 48 pair 
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Table 2. Basic information for a particularly striking example of 
a Milky- Way Magellanic Cloud Analogue (GAMA-MMA1) with 
a near-by bright spiral companion. This mimics a lot of the most 
recognisable characteristics of the MW-SMC-LMC and M31 com- 
plex. 

systems). For each pair system we calculate the effective r- 
band centre-of- light in RA, Dec and redshift and define this 
as the centre of the pair system. We now search for all MCAs 
that are within a r pro j < l,000kpc and v sap < 500 km s 
separation to L* pair systems. 11/47 MCAs are found in 
close proximity to L* pair systems, while 27/47 MCAs are 
found within the same spatial separation of any L* galaxy. 
At a maximum the effective comoving volume searched over 
for all L* systems is 1.3 x 10 4 Mpc 3 , which is ~7% of the 
available GAMA volume within these redshift limits. For the 
L* pairs the maximum volume searched is 2.1 x 10 3 Mpc 3 , 
which is ~1% of the available GAMA volume. Consequently 
we find that 57% of MCAs are found within 7% of the vol- 
ume when searching around L* systems, and 23% of MCAs 
are found within 1% of the volume when searching around 
the centres of L* pairs. There is a clear tendency for MCAs 
to be associated with L* systems in general (rather than 
just random distributed throughout the Universe), and an 
even stronger association is seen between MCAs and binary 
pairs of L* galaxies, like the MW and M31. 

The best MMA (constituting a spiral brightest-galaxy 
and star forming companions, known as GAMA-MMA1) 
also has a near-by companion spiral galaxy. A multi- 
colour image based on SDSS photometric data is shown 
in Figure [2] Table [2] contains key information on the sys- 
tem. The companion spiral is 190 Kpc away in projec- 
tion (MW and M31 are separated by 800 Kpc in real 
space) and -400 kms _1 separated in velocity (M31 is mov- 
ing at — 122 km s _1 radially and ~ 100kms _1 tangentially 
with respect to the MW). These two spirals have M r — 
-21.43/ - 21.31, which is very similar to MW / M31 (M r = 
—21.47/ - 21.17). IFU and Parkes data have recently been 
obtained for this newly identified analogue to the LG sys- 
tem, and we are actively seeking followup on a sample of the 
MMAs discussed in this work in order to better quantify the 
occupation statistics in halos similar to our own. 

5 CONCLUSIONS 

The major findings of this work are summarised below: 

• Analysing all galaxies within 0.3 dex of the stellar mass 
of the MW, there is a 11.9% (U.2%-12.8%) chance that it 
will have a close companion (within a projected separation 
of 70 kpc and radial separation of 400 km/s) that is at least 
as massive as the LMC. This is consistent with analyses by 
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Figure 2. SDSS image of the best LG analogue found in GAM A. 
The spiral galaxy in the bottom right with the two indicated 
brighter companions is the MMA (GAMA-MMA1), where all 
three galaxies are late-type and star-forming, making this system 
similar to the MW-Magellanic Cloud system in the most funda- 
mental ways. An associated nearby spiral galaxy is also shown. 



iBovlan-Kolchin et alJ ( 201 ll) of the Millenium II simulation 
and by James fc Ivory ( 201lh of Ha imaging around lumi- 
nous spiral galaxies. 

• Limiting the sample to those galaxies where the SMC 
would be observable we find 3.4% (2.7%-4.5%) of galaxies 
have two companions at least as massive as the SMC. 

• Only two full analogues to the MW-LMC-SMC system 
were found in GAM A, suggesting such a combination of late- 
type, close star-forming galaxies is quite rare: in GAMA 
only 0.4% (0.3%-l.l%) of MW mass galaxies have such a 
system (a 2.7a event). In terms of space density, we find 
1.1 x 10 -5 Mpc~ 3 full analogues in GAMA (in a volume of 
1.8 x 10 5 Mpc 3 ). The best example found shares many qual- 
itative characteristics with the MW system. The brightest 
pair galaxy has spiral features, as does the bigger minor 
companion. The minor companions are ~40 kpc in projected 
separation, so not in a close binary formation like the SMC 
and LMC. 

• Selecting systems that are close binaries like the SMC- 
LMC pair (MCAs), we find that they are preferentially lo- 
cated in close proximity (or within) systems that have a 
similar total flux to the LG (SM r = -22.5 ± 0.5 mag). 

• Subdividing the preferential group type into those with 
large and small magnitude gaps, we find MCAs are more 
spatially associated with groups that have a small magni- 
tude gap. This suggests that a quiet recent merger history 
improves the likelihood of the Magellanic Clouds being vis- 
ible in the LG. The best MMA analogue found in GAMA 
also has a close L* spiral companion galaxy. 
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